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Neumann 问题 热平衡 积分 法 细 化 解 的 收敛 性 * 
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摘 要: 本 文 应 用 热平衡 积分 法 及 其 细 化 方法 求 得 了 Neumann 问题 的 分 段 线性 近似 解 。 在 参数 S 充分 大 
的 条 件 下 ， 应 用 伽 玛 函 数 及 不 完全 伽 玛 函数 的 性 质 ， 证 明了 Neumann 问题 的 热平衡 积分 细 化 解 
的 收敛 性 ， 并 求 得 各 相 的 收敛 速度 分 别 为 DO(n- 0)、O(n-22)。 本 文 推广 了 热平衡 积分 细 化 解 收 
敛 性 分 析 中 的 一 些 已 知 结论 。 
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热平衡 积分 法 (Heat Balance Integral Method, HBIM) Æ Goodman!) 于 1958 年 提出 的 求解 
导热 问题 的 一 种 近似 方法 。 大 量 实例 表明 ， 选 取 适 当 的 近似 函数 就 能 得 到 所 求 问题 的 较 好 的 
近似 解 。 然 而 ，LangfordD] 举例 说 明了 HBIM 中 近似 函数 的 选取 无 统一 的 办 法 ， 并 指出 使 用 
高 次 多 项 式 不 一 定 能 提高 近似 解 精度 。 为 解决 这 一 问题 ，Noblegl 提出 用 剂 分 积分 区 间 与 分 段 
低 次 多 项 式 相 结合 的 办 法 ( 即 细 化 、Refinement)。Bell 在 文献 4,5] 最 先 应 用 了 该 方法 ， 其 数值 
结果 证 实 了 HBIM 细 化 的 有 效 性 。1985 年 ， 文 献 [6] 理论 上 证 明了 半 无 限 纯 导 热 问题 HBIM 细 
化 解 的 收敛 性 。 但 仅 限于 纯 导 热 问题 ， 未 讨论 相 变 问题 。2005 年 ，Mosally 等 人 结合 数值 实验 
给 出 了 单 相 融化 问题 的 HBIM 细 化 解 在 Stefan 数 8 = 1 时 的 收敛 速度 由， 但 其 结果 依赖 于 数 
值 实验 ， 故 不 能 作为 严格 的 理论 依据 。2008 年， 文献 [8] 证 明了 单 相 融化 问题 的 HBIM 细 化 解 
É Stefan & GB > 1 时 收敛 于 精确 解 ， 给 出 了 收敛 速度 。 

工程 实践 中 导热 问题 大 多 是 两 相 或 多 相 的 ， 一 般 难 以 求 其 解析 解 ， 人 们 通常 用 近似 方法 求 
解 此 类 问题 。 本 文 从 理论 上 考察 将 HBIM 和 细 化 方法 应 用 于 这 类 问题 的 合理 性 ， 为 便于 比较 ， 
文中 以 解析 解 存在 的 Neumann 问题 四 作为 两 相 问题 的 模型 ， 结 合 文献 [6,7] 中 的 方法 ， 给 出 了 
该 问题 的 热平衡 积分 细 化 解 的 收敛 条 件 和 收敛 速度 。 


2 Neumann 问题 的 热平衡 积分 模型 
Neumann 问题 区 的 数学 描述 如 下 


一 - 一 < 

Bt = gra 0<z< X(t), (1) 
OT» O°T2 

—_—_— = aes =p < 

BE = 02 Ba X(t) <x < +00, (2) 
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T,(0,t) = Tos t> 0, (3) 
im To(x,t)=T,, t>0, (4) 
Ty(X,t) = To(X,t) = Tf, (5) 
OT, OT» a dX 加 
ha koa pl r=X (6) 
(0, X] 为 第 一 相 ，[X, +00] 为 第 二 相 ， 
D1C1 P2C2 á 
pi ci (i = 12) 分 别 为 两 相 的 热传导 系数 、 密 度 、 热 容 。aa, as 为 两 相 热 扩散 系数 。 记 ai/a2 = 
arz, ka/kı = kz1， 已 知 问题 (1)-(6) 的 解析 解 四 为 
x 
T =T, +A ert i) (7) 
T, = T.~ B[1—erf( z7 )| =T- Bete 5 ma) (8) 
X = 2% Vout, (9) 
其 中 erf(n) 称 为 误差 函数 ， 其 表达 式 为 
PEME 
erf (7) = = | e* dé, 
erfc (n) = 1 一 erf (n) 称 为 误差 函数 的 余 函 数 。 
-Dan po ni 
erf(yY)’ © erfe(yyar2) 
7 为 如 下 超越 方程 的 根 
e kavan Tyee" lyya (10) 
erf(y) (Ty ~Ts)erfc(y/oi2) ca — Ts) 
由 解析 解 易 知 i 
—2=0, t>0. (11) 


为 简化 计算 ， 在 问题 (1)-(6) FR 
Ts=-1, Ty=0, To =1, 


FA SCAR [6,7] I, Ty = OW, Ty, To AER MAA MBIT, AAA 
HIT, To» IEA AF ANA TAR REA). HT, = 0 总 能 做 到 ， 若 Try £0, (RBH 
HT =T —T;, Th =To —T;, R (5) KALA 


Ty (X, t) = T,(X,t) - Tr =0, Ty(X, t) = To(X,t) — T; =0, 


Aid: S= l/cio 
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经 上 述 简化 后 ， 根 据 HBIM KEN Xf Neumann 问题 建立 热平衡 积分 模型 ， 即 应 用 热平衡 
积分 法 对 Neumann 问题 进行 近似 求解 。 对 半 无 限 问题 应 用 HBIM 需 引 入 摄 动 深度 四， 即 在 上 时 
刻 ， 假 设 固 定 边界 z = 0 处 的 温度 扰动 只 影响 到 有 限 的 深度 6(t) 处， 超过 这 一 深度 处 物体 的 温 
度 仍 为 初始 温度 7 ， 也 就 是 热 交 换 不 再 发 生 。 这 就 是 对 温度 扰动 的 影响 在 瞬间 达到 无 穷 远 处 的 
近似 。 并 从 Neumann 问题 (1)-(6) 解析 解 的 角度 分 析 四 有 


< =b, (12) 


5b 为 与 时 间 无 关 的 常数 且 b > 1。 引 入 摄 动 深度 56( 后 ， 则 在 上 + 时 刻 时 只 需 考虑 区 闻 [0, 8t) 内 的 
温度 分 布 ， 于 是 对 第 二 相 控制 方程 的 求解 区 间 和 边界 条 件 有 如 下 近似 


X<zr<ð, (13) 
T(t =T, =1, t>0. (14) 


将 (13)、(14) 代入 (2); (4) 得 到 问题 (1)-(6) 的 HBIM 求解 模型 。 


3 ” 细 化 模型 和 分 段 线性 热平衡 积分 解 


KT, Te 的 空间 变量 定义 域 [0,X], 区 ,6] 各 自 等 分 成 有 份 ， 并 记 


iX i(6—X) _ 
一 一- 


X,=—, Y=X+ pe ee i=0,---,n. 
n n 


TER X b-1)X 
Xini -Xi = >, preke dh 
n n 


WU}, US DIA Ti, Tz 在 节点 i 处 的 近似 值 ， 由 边界 条 件 有 
U? = —1, Ur=UV9=0, Ur?=1. 

在 [Xi, Xi] Lie ERATED T (x,t) 的 近似 函数 ， 记 为 
n(Uit’ — Ui)(z — Xi) 
x f 


Uı = Uj + z € [Xi, Xiz], 1=0,--- ‘n=l, (15) 
在 [Y;, Yip] LEARE BE To(x, t) 的 近似 函数 


n(U3** — U3)(z ~ ¥i) 


U, = U% + GX ， ©€(%, Yui], i=0,..……,n—l1, (16) 
由 (15)、(16) 得 
aU, — n(Ui*1 — Ui 
ad we 2: zE [Xi Xim) 1=0,---,n—-1, (17) 
i+l _ rr 
OU2 _ n(U3 ey zE [Yn Yn) i=0,---,n—1, (18) 


ðr — (b-1)X 


vit _ui=Ai, U$ -Ui = A$, 
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在 [Xi, Xim], [Yi Yiyi] (i = 0,… ,n 一 2) 上 分 别 对 (1) 式 和 (13) 式 进行 热平衡 积分 ， 根 据 莱 布 
尼 茨 法 则 并 用 U RET, Uo RET, BHA 


dX 2n2a (Ai—Ait!l) . 
= y : 22 Qe ese —2 1 
dX _ 2n? a2 (A2 一 Ai ) L 
ET 
同 理 ， 在 [Xn_1, Xn], Yn, Yn] 上 对 各 自 的 方程 进行 积分 得 
dx = 2n?a [(b 一 DA! = ko, A9] (21) 
dt — (b—1)[2nS + (2n—1)A?“}]’ 
2 
xi 2n a2 (22) 


a (b—1)(2nb—b+1)’ 


其 中 在 (22) RFE SIE PER Yn = (t) 处 的 温度 梯度 

OU po 

Bz È t) =0, 
这 是 基于 超过 摄 动 深度 处 不 再 发 生 热 交换 的 假设 ， 也 就 是 温度 梯度 为 零 。 同 时 ， 作 为 近似 解 ， 
这 也 和 解析 解 的 性 质 (11) 式 相 符 。 令 


= 


5 n+1/2, (23) 


H (22). (23) 得 
n 
X(t) =2- Vent ont vart, 
记 
Y = (24) 


VF I 
又 由 (19)、(22)、(23) 联 立 ，(20)、(22)、(23) 联 立 可 得 


i . i+1 一 see a 
Ai 和 Ay, i=0,---,n-2, (25) 
t @-1/2 a 
i b-1 . Ait oe | ee 到 
入 = A i=0,.,n-2. (26) 
另外 有 如 下 关系 

n-1 n-1 
$ Af=up-ups=1, So Ab =up -U =1, (27) 
k=0 k=0 


PRUE (25). (26). (27) 可 解 得 含 参 数 5 的 AY, AS (i = 0,.… ,n 一 1)， 注 意 到 (21)、(22) 有 


an[(b — 1)A? ’ — k21A9] 加 Q2 (28) 
QnS + (2n — 1)A7-! 2nb—b+1’ 
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FER GH SR OK AR, ASTRA (26) 可 求 得 6， 从 而 求 得 的 表达 式 和 Ai, Ab (i = 
0,… ,n 一 1)， 所 以 解 得 Ui, Ui (G=1,---,n-1) 


n—l n—l 
Ui=U?r—D At=— A‘, U= Up- -1 (29) 
k=i =i k=i 


4 ”收敛 性 分 析 


问题 (1)-(6) 的 热平衡 积分 细 化 解 的 收敛 性 分 析 ， 需 要 导出 Ai, Ui 和 人 3, US 的 具体 表达 
式 。 由 (25)、(26)、(27)、(29) 并 应 用 伽 玛 函 数 的 性 质 F(z + 1) = zT(z) 可 得 


2 1 
An} = = (30) 
P(n + 1) 2 下 


(ntn—1/2) -iD 
Ai = er) 一， i=0,.. ,nl (31) 

> (n+n—1/2)* 

er T(n+k+1) 


mA (n+n-1/2) 
P(n+k+1) 


Uf SR isn], (32) 
nt+n~1/2)* 
res 


1 
ae ca ee (33) 


(bn/(b—1)-1/2)* 
T(k+1 


k=0 


(bn/(b-1)-1/2)"7* 
AS TO T(n 


n~-1 ? 
(6n/(6—1)—1/2)* 
和 T(k+1) 


(34) 


全 nor ae 1a 


Ui = 1— + ， i=0,; n=l. (35) 
2 utie Dn =a 


又 由 (10) 式 易 知 ， 对 充分 大 的 S， 则 ?7 很 小 ， 于 是 对 这 样 的 7Y， 误 差 函数 erf(7) 有 渐 近 
式 [19] 


erf (7) = -= -e77 - [y+0(7)), 


VT 
将 其 代入 (10) 式 并 忽略 Y 的 二 阶 以 上 的 无 穷 小 量 解 得 
1 
1 (36) 


当 5 充 分 大 时 应 用 (36) 式 可 得 如 下 引 理 。 
引 理 当 5 充 分 大 时 ， 对 (24) 式 中 的 YY， 有 lim y = Y， 其 收敛 速度 为 On )。 


508 I E 数 学 学 R 第 27 卷 


证 明 注意 到 (30) 式 中 分 母 中 的 每 一 项 


(n+n-— 1/2) = bes ate) 
Na+) Trit 2 k=0,- n 1, 
故 有 
(n+ ~ 1/2)* 
roy ar ey 
PRA AT + <4, HH (28) RA 
Q2 ag ab 


Inb ~ Inb—b+1 ~ 2n2S” 


从 而 有 
> 
Q12 
即 5 发 散 于 正 无 穷 。 因 此 当 n 一 +oo 时 ， 有 5b/(b 一 1) 一 1。 所 以 当 n 一 +co 时 ， 由 (34) 式 知 
(n=1/2)"—! 
A}? — T(n) 


n—l i 
(n—1/2)* 

X Gu 

再 根据 Bell 和 Abbas X PHRAS, 4n 一 +oo 时 ， 


(n — 1/2)"-1/2 E 
? 三 -站 
"Sl (n-1/2)* T 


T(n) D TFI 


于 是 得 
(n—1/2)"7? 5 
Ag ,— rín) 5 ae (37) 
ee 
又 由 5b 发 散 于 无 穷 大 ， 故 n 一 +00 时 
7 1 一 1/2 > n-1 1 1 i 2 
T= ane( SP - Rb 4 =a b*, (38) 
H (37). (38) 代入 (28) 式 中 并 略 去 无 穷 小 量 得 
a 
Al = Am (39) 


其 中 c= 2. k3 -ai M (38) R4 n EEK, n> Sn2， 也 即 
n 
< /3 < Vn, 


(n+n— 1/2)* 1 k 
( 1) T( k 1) <(1 I > 区 1, 
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所 以 当 5 充 分 大 时 ， 这 里 指使 得 
时 ， 有 
n-1 n-1 
(n +7 — 1/2): ey b EE lk 
T+) Sarry < (+t) Se Dy, 


所 以 由 上 式 结合 (30)、(39) 式 得 


对 充分 大 的 S 有 


(y5) 


由 (40) 式 有 7 < n25， 因 此 得 n= nmn2S。 于 是 有 


lim y= 
cares? 


k n ne Pee 
nut, J2n + 2n—1 ~ V25 TN 
BAY =7+O(n""). 证 毕 
为 考察 边界 x = 0 处 的 热流 和 节点 处 温度 的 收敛 情况 ， 引 入 不 完全 伽 玛 函数 D 


,Z) 一 er-ttar-ldt 
g(a, z) f e 
由 分 部 积分 得 
g(a,2) = (a — 1)g(a — 1, z) 一 ez 1. (41) 
对 边界 z = 0 处 的 近似 热流 ki ZA 和 近似 温度 Uj A Fe 
定理 1 4SRAAN 


aU, 
noo Gr T! 


且 二 者 的 收敛 速度 均 为 O(n-!)。 
证 明 由 (17) 和 (31) 式 得 


1 r lz=0 ~ 2y* Vat Viait k 
i i rin +n) E E 


z=0, lim Uj =T}, 
n~— +00 


又 由 引 理 和 Mossally 等 人 文中 的 结论 中 有 


(nt+n—1/2)7"7? V2 


oo (a -1/2  Vrerf(y)’ 
P(n +n) 2, O 


lim 
n+ 
并 知 其 收 钱 速度 为 O(n-!)。 所 以 证 得 


imk ðU: = ky =f. OT; 
nto! Oz le=0 /moyterf(y) | Ox le=0’ 
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RA O(n!) 的 收敛 速度 。 再 根据 递 推 关系 (41) 可 将 (32) 式 化 为 


g(nn+n-1/2) _ g(nty—i,ntn—1/2) 
T(n) n+n—i 


Uis g(nn+n—1/2) _ g(ntn,n+n—1/2) ’ 
T(n) T'(nt+n) 


式 (42) 右边 分 子 第 一 项 和 分 母 的 收敛 性 由 引 理 和 Mosally 等 人 中 的 结论 可 知 


St 一 -31i+eto+of 二 )]， 


T(n) 2 vi 
g(in+nn+n-1/2)_1 1 
Pren) =3[1+0( a) I: 
对 于 (42) 式 分 子 第 二 项 的 收敛 性 ， 注 意 到 对 固定 的 点 (x,t) 有 
eee 
n 


ix (i+1)X 
Va nyt 
所 以 当 一 上 occ 时 ，i ~ EV2(n 十 9) 从 而 得 


g(n+7—int+n— 1/2) = 5 [ite +0( =—)|. 


Ynt+n-t 


T(n+7-2) 
所 以 由 (43). (44), (45) 三 式 证 得 
E ts Gea Tare) Se 
ee ae et 
并 易 知 其 收敛 速度 为 O(n-1/2)， 也 就 是 O(n-!1)。 
MF US 的 收敛 性 有 如 下 定理 。 
定理 2 ” 当 5 充 分 大 时 有 


lim Ui= Ti 
A 2 2: 


其 收敛 速度 均 为 O(n-1/2)。 
证 明 由 递 推 关 系 (41) 可 将 (35) 式 转换 成 如 下 形式 


i= g(n~i,bn/(b—1)—1/2) 
Pai 
2 1 一 ECD2I3 ， 
rín) 
又 由 Mosally FAKAT, WEEKK mA 


aa = [1 + erf( is ) + O(m-¥?)], 


应 用 (47) 式 和 


1267 = 7 = ns, 


(42) 


(43) 


(44) 


(45) 


证 毕 


(46) 


(47) 
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的 关系 得 出 
sui O—D =U) 1 4 gt ODaW=*) 4 Og] 
= 5[ttert( vt) + O(n-¥)], 


Xf Te E NRH IA FE AY xa (a, t) A 


x 4 OK cg cx OEL, 


同样 取 é& = Want’ 故 依 照 定 理 1 中 的 方法 知 ， 当 ?mm 一 十 co 时 


i~ adc ee Ga 
b—1 
所 以 由 (47) 式 和 4b 的 性 质 
g(n — i, bn/(b— 1) — 1/2) _ 1 sai) 
to )] 
因此 由 (46) 式 得 
| ; ， iret) efc(sg) ,, 
gee) Taig) yaa) 
且 其 收敛 速度 为 n-1/2。 证 毕 
5 it 


本 文 应 用 HBIM 的 细 化 求 得 Neumann 问题 的 热平衡 积分 细 化 解 ， 证 明了 该 近似 解 在 参 
数 5 充分 大 时 收敛 到 解析 解 ， 并 求 得 各 相 的 收敛 速度 。 和 两 相 问 题 类 似 ， 本 文 的 细 化 模型 和 求 
解 方法 可 用 于 多 相传 热 问题 ， 文 中 收敛 性 结果 为 HBIM 细 化 法 在 多 相 问 题 中 的 应 用 提供 了 理论 
依据 。 此 外 ， 本 文 推广 了 文献 [6,7] 中 的 结论 。 事 实 上 ， 可 以 把 Neumann 问题 的 第 一 相 看 作 单 
相 变 问题 ， 第 二 相 看 作 纯 导热 问题 。 纯 导热 问题 则 可 将 摄 动 深度 视 为 相 变 面 ， 故 可 看 作 相 变 问 
题 来 处 理 。 另 外 ， 分 析 文 献 [6,7] 和 本 文 的 结论 ， 发 现 融 化 参数 y 的 收敛 性 质 决 定 其 它 量 的 收敛 
EH, KER y 的 收敛 性 与 温度 场 及 热流 的 收敛 性 之 间 的 内 在 联系 应 该 在 HBIM 的 研究 中 引起 
重视 。 


参考 文献 : 


[1] Goodman T R. The heat-balance integral and its application to problems involving a change of phase[Jj. 
Tran ASMEJ Heat Transfer, 1958, 80: 335-342 

[2] Langford D. The heat balance integral method[J]. Int J Heat Mass Transfer, 1973, 16: 2424-2428 

[3] Noble B. Heat balance methods in melting problems[C]// Ockendon J R, Hodgkins W R (eds.), Moving 
Boundary Problems in Heat Flow and Diffusion, Oxford Conference, March 1974, Oxford: Clarendon Press, 
1975 


512 I E 数 学 学 报 第 27 卷 


[4] Bell G E. A refinement of the heat balance integral method applied to a melting problem[J]. Int J Heat 


Mass Transfer, 1978, 21: 1357-1362 
[5] Bell G E. Solidification of a liquid about a cylindrical pipe[J]. Int J Heat Mass Transfer, 1979, 22: 1681-1686 
[6] Bell G E, Abbas S K. Convergence properties of the heat balance integral method[J]. Numer Heat Transfer, 


1985, 8: 373-382 
[7] Mosally F, Wood A S, AL-Fhaid A. On the convergence of the heat balance integral method|[J]. Applied 


Mathematical Modelling, 2005, 29: 903-912 
[8] 徐 湘 田 ， 令 锋 . 单 相 融 化 问题 热平衡 积分 细 化 解 的 收敛 性 [了 ]. 内 蒙古 大 学 学 报 (自然 科学 版 )，2009, 40(2): 6-9 
Xu X T, Ling F. Convergence properties of the refined heat balance integral solution for the set of one-phase 
melting broblem[J]. Journal of Inner Mongolia University, 2009, 40(2): 6-9 
[9] Lunadini V J. Heat Transfer in Cold Climate[M]. New York: Van Nostrand Reinhold Company, 1981 
[10] Wood A S. A new look at the heat balance integral method[J]. Applied Mathematical Modelling, 2001, 25: 


815-824 
[11] Abramowitz M, Stegun I A. Handbook of Mathematical Functions[M]. New York: Dover Press, 1965 


Convergence Properties of the Refined Heat Balance Integral 
Solution for the Neumann Problem 


XU Xiang-tian’?, LING Feng? 


(1- College of science, Inner Mongolia University of Technology, Hohhot 010051; 
2- School of Computer Science, Zhaoqing University, Zhaoqing 526061) 


Abstract: A piecewise linear approximate solution for the Neumann problem is presented by using 
the heat balance integral method and its refined method. It is proved that the approximate solution 
converges formally to the exact solution at large values of the parameter S and the convergence rate 
in each phase is O(n~') and O(n! 2), respectively, by using the properties of the Gamma function 
and incomplete Gamma function. The methods used in this paper can be applied to two-phase and 
multi-phase heat conduction problems. The present work generalizes the former results in convergence 
analysis to the refined heat balance integral solution. 
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